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ABSTRACT 

We present spherically symmetric and plane parallel photoionization models 
of NGC 346, an H II region in the Small Magellanic Cloud. The models are based 
on Cloudy and on the observations of Peimbert, Peimbert, & Ruiz (2000). We 
find that approximately 45% of the H ionization photons escape from the Hll 
region providing an important ionizing source for the low density interstellar 
medium of the SMC. The predicted J(4363)//(5007) value is smaller than that 
observed, probably implying that there is an additional source of energy not taken 
into account by the models. From the ionization structure of the best model and 
the observed line intensities we determine the abundances of N, Ne, S, Ar, and 
Fe relative to O. 

Subject headings: galaxies: abundances — galaxies: individual (SMC) — galaxies: 
ISM — H II regions — ISM: abundances 



Introduction 



NGC 346 is the most luminous Hll region of the SMC, its Ha luminosity places it 
on the boundary between normal and giant extragalactic Hll regions. We have produced 
photoionization models of this object for the following reasons: a) to find if NGC 346 is 
density bounded or ionization bounded, b) to determine its chemical composition, and c) 
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to study the energy budget of the region. The photoionization models have been computed 
with Cloudy 94 (Ferland 1996; Ferland et al. 1998). 

The photoionization models computed in this paper present two advantages relative to 
similar models of extragalactic H II regions available in the literature: a) there are excellent 
line intensity observations available for this Hn region (Peimbert, et al. 2000, hereinafter 
Paper I) ; b) the spectral types of most of the ionizing stars are known which permits us to 
adopt a representative ionizing flux. Most photoionization models in the literature are based 
on three assumptions: an IMF, an analytical fit to the adopted IMF, and a star formation 
history, these three assumptions introduce significant uncertainties in the ionizing flux (e.g., 
Cerviho, Luridiana, & Castander 2000). 

2. Stellar Content and Photoinization Spectrum 

We estimate directly the photoionization spectrum that will be used for the models 
based on the spectral classification of the blue stars presented by Massey, Parker, & Garmany 
(1989). 

The ionized region of NGC 346 has a diameter of about 420" and from the figures and 
coordinates presented by Massey et al. (1989) and Ye, Turtle, & Kennicutt (1991) we find 
that there are 58 blue stars with a visual apparent magnitude V < 15.6 inside the ionized 
boundaries of NGC 346. The blue stars fainter than V = 15.6 are expected to be of spectral 
type 09 or later and do not contribute appreciably to the photoionization spectrum. 

Of the 58 stars considered by us 33 have pre-assigned spectral types and 25 do not. We 
have counted binary systems as single stars with the exception of HD 5980 where we consider 
independently the two main components (the system might have a third component). 

We have assigned spectral types for the other 25 objects based on their V magnitudes 
and assuming that they are on the main sequence, these are the objects in parentheses 
presented in Table 1. 

To derive the ionizing flux for all the stars in Table 1, with the exception of HD 5980, 
we have adopted the parameters presented by Vacca, Garmany, & Shull (1996). 

For HD 5980 we have taken the stellar parameters derived by Schweickhardt & Schmutz(1999) 
where they have assumed that the luminosity of each component has remained constant dur- 
ing the 1987-1999 period, but varying the temperature and radius of component A. For the 
temperature of component A we have adopted the value it had at the end of 1990, the time 
of the observations that we are trying to match with the models. At this temperature the 
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stellar radius had the value of 21.4 R . 

The total ionizing flux is 40.08 x 10 39 erg s _1 (see Table 1). The ionizing fluxes of the 
objects without spectral types are presented in parentheses in Table 1; their combined flux 
amounts to 8.80 x 10 39 erg s -1 , 21.9% of the total value, which is a significant but relatively 
small fraction. 

We estimate errors of 0.2 dex, 0.2 dex and 0.08 dex for the flux of HD 5980, the total 
flux of the stars in parentheses, and the total flux of all the other stars with known spectral 
types, respectively; consequently we have adopted an error of 0.08 dex for the total ionizing 
flux. 

From the ionizing flux and the effective temperature for each star presented in Table 1 
we have built three different photoionization spectra by adopting : a) blackbodies for all the 
stars, b) a set of stellar models for = Z Q , and c) a set of stellar models for Z* = 0.2 Z & . 
These three ionization spectra were chosen to estimate how sensitive are the properties of 
the models to changes in the adopted radiation field and consequently how general or robust 
are the results obtained from the models. The low metallicity set is based on the gaseous 
abundances for NGC 346 derived in Paper I and assumes that 20% of the heavy elements 
are trapped by dust inside the Hn region (Esteban et al. 1998). 

For the models with Z* = Z@, we chose: Mihalas (1972) atmospheres for the 03V star 
and for 5980 A, Kurucz (1991) atmospheres for the stars of Luminosity Class I, and Schaerer 
& de Koter (1997) for the rest of the stars in Table 1. The models with Z* = 0.2 Z Q are 
obtained using the library of stellar spectra built by Lejeune, Cuisinier, & Buser (1997); 
this library is based on three sets of atmosphere models and covers a wide range of physical 
stellar parameters. The three different ionizing continuum spectra used for the models are 
illustrated in Figure 1. The model sets are named as followed; series of models B, for 
blackbody spectra, series of models S for stellar models with Z* = Z & and series of models 
L for stellar spectra from Lejeune, Cuisinier, & Buser (1997). 

Based on the models just described and assuming that all the ionizing photons are 
trapped by the nebula, the total ionizing flux (that is the same for all the models) implies 
an Ha flux, L(Ha) in erg s" 1 , of 39.216 dex for the model with Z* = Z Q , 39.237 dex for 
the model with Z* = 0.2 Z & , and 39.25 dex assuming a blackbody spectrum. The small 
differences in the Ha flux are due to the difference in the number of ionizing photons, which 
for a given ionizing flux does depend on the shape of the spectrum. 

On the other hand, from the observed flux of 1.55 x 10 -9 erg s _1 cm -2 (Kennicutt & 
Hodge 1996), an interstellar absorption correction, C(Ha), of 0.10 dex (Paper I), and a 
distance to the SMC of 64 kpc (Reid 1999), we obtain a total emitted Ha flux of 38.98 ± 
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0.06 dex (erg s 1 ). Where the error is due to the combination of 0.04 dex, 0.02 dex, 0.01 dex 
errors in the observed flux, the adopted distance and C(Ha), respectively. 

The difference between the Ha flux predicted from the ionizing stellar radiation and 
the Ha flux derived from observations amounts to 0.24 ± 0.10 dex, 0.26 ± 0.10 dex and 
0.27 ± 0.10 dex for the models with Z* = Z Q , Z* = 0.2 Z Q and for a blackbody spectrum 
respectively. These values imply that about 45 ± 15% of the ionizing photons escape from 
NGC 346, and that the region is density bounded. Moreover, this result is in excellent 
agreement with the models to ionize the diffuse interstellar medium proposed by Zurita, 
Rozas, & Beckman (2000), and with the results of Oey & Kennicutt (1997), who find that 
many Hll regions in the Large Magallanic Cloud are density bounded. This problem is 
discussed further in section 3.6. 

The result that 45% of the photons escape is independent of the electron density dis- 
tribution, but it does not necessarily imply that the Hll region is density bounded in all 
directions, it is also posible to have a region which is density bounded in some directions and 
ionization bounded in others, these types of models require density distributions that are 
not radially symmetric, the simplest models that study this effect are those with a "covering 
factor", where in some directions all the photons escape and in others all the photons are 
trapped. We have assumed for the spherical models a constant density distribution or a 
density distribution that is spherically symmetric, therefore by density bounded we mean 
that photons escape in all directions. It is beyond the scope of this paper to model other 
types of density distributions. 



2.1. Initial Mass Function 

It is important to compare our ionizing flux distribution with that derived from a 
Salpeter IMF, for this purpose we will estimate the IMF of the massive stars responsible 
of most of the ionization. 

In Figure 2 we present the IMF derived from Table 1 under the following assumptions: 
a) the main sequence stellar masses have been adopted as the initial masses, b) the mass of 
stars that are not in the main sequence has been adopted as their main sequence mass, c) for 
each star instead of an exact mass value we have adopted a uniform continuous distribution 
given by the central value ± 0.02 dex, d) we have fixed the lower mass limit of the last 
bin at 51 M , which implies that this bin includes only half a star, e) the width of the bins 
corresponds to 0.08 dex, f) to derive the slope of the IMF we only considered the central four 
bins, we did not include the last bin because the IMF is truncated at the high mass end which 
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is inside this bin. We did not include the bin at the low mass end due to incompleteness, 
moreover the fraction of the ionizing flux due to the stars in this bin is negligible. 

The fit to the central four bins defined in Figure 2 is given by 

log[N/A log m] = 7.675 - (3.71 ± 0.4) log m (1) 

this equation and point d) of the previous paragraph implies that the upper mass cutoff 
is equal to 54.1 M . The slope given by equation (1) is considerably steeper than -2.35, the 
Salpeter slope, and indicates that for NGC 346 the use of the Salpeter slope to represent the 
ionizing flux is not adequate. 

There are two other determinations of the IMF for the SMC. Massey et al. (1989) obtain 
-2.9 ± 0.3 in the 9 M to 85 M mass range for NGC 346 and Humphreys and McElroy 
(1984) obtain -3.1 in the 15 M to 100 M mass range for the SMC. From the data of 
Humphreys and McElroy we obtain a slope of -3.7 for the 25 M to 100 M mass range, in 
excellent agreement with the slope derived by us. 

Our IMF for NGC 346 implies that for extragalactic H II regions with an Ha flux similar 
or smaller than that of NGC 346 it is difficult to determine M up , and that the slope, at least 
for masses higher than 24 M , might be steeper than the slope given by a Salpeter IMF. 

To determine the ionizing flux for an extragalactic Hll region based on an analytical 
function we need to know: a) the slope of the IMF, b) the M up , c) the fraction of ionizing 
photons that are trapped by the H II region, and d) the star formation history of the burst 
that produced the Hll region. As discussed above these four problems indicate that for 
NGC 346 it is considerably better to use the ionizing flux determined from the individual 
stars than to derive the ionizing flux based on an analytical approximation to the IMF. 

3. Photoionization Models 

All the models discussed in this paper were computed with Cloudy. For all models, 
with the exception of those with enhanced and depleted abundances (models L.5 and L.6), 
we adopted the following chemical composition (given in log A(X)/iV(H) + 12): He = 10.90, 
C = 7.39, N = 6.51, O = 8.11, Ne = 7.30, Mg = 5.98, Si = 6.37, S = 6.59, CI = 4.80, Ar = 
5.82, and Fe = 5.58. Where the He, N, O, Ne, S and Ar are the gaseous abundances derived 
in Paper I, C/O is the value derived by Garnett et al. (1995) for the SMC, and Mg/O, Si/O 
and Fe/O come from the gaseous abundances derived by Esteban et al. (1998) for the Orion 
nebula. 
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As mentioned in section 2 we used three different sets of ionizing ffuxes: atmospheres 
with = 0.2 Z Q (models L), atmospheres with = Z & (models S), and blackbody 
radiation fields (models B). The reason to use these three sets of ionizing fluxes is to explore 
how general are the conclusions based on the L models. In Figure 1 we show the three 
integrated spectra used as the input ionization for the models. 

To compare the observed values with the models computed with Cloudy we used the 
radial integrations of the models that include the center and compared them with the average 
intensities of Regions 3 and 13 which are close to the center. We took the average of the 
two regions to reduce the observational errors and to obtain a more representative value of 
a line of sight towards the center of the object. All the models have an inner radius, ro, of 
10 16 cm, with the exception of models L.7 and B.3. 

We also present in Tables 2 and 3 the line intensities of region A (that includes the 
average intensities of 13 regions) because the errors are smaller than those of the average 
intensities of regions 3 and 13 and to show that the differences between the two sets of line 
intensities are relatively small. For similar reasons we also present T e (rad)K and T e (vol)K the 
radial and the volume averaged temperatures weighted by N e N{0 ++ ) obtained from the mod- 
els, and T e ([Olll])K the volume averaged temperature obtained from the J(4363)//(5007) 
ratio. To compare the line intensities with those of regions 3 and 13, all the line intensi- 
ties have been radial averaged. The main reason for selecting the line ratios presented in 
Tables 2, 3 and 4 is that they are very sensitive to changes in the input parameters for the 
models, thus allowing us to obtain a robust model for the observed region. In particular the 
[O n]/[0 ill] ratios are very sensitive to the degree of ionization, and the [O il], [S n], and [O i] 
to H/9 line ratios give us information on the behaviour of the low degree of ionization regions 
of the nebula. All these line ratios are also sensitive to the density distribution adopted. 

In sections 3.1 to 3.5 we discuss density bounded models in all directions in agreement 
with the results presented in section 2. In section 3.6 we discuss two ionization bounded 
models which provide independent support for the view that NGC 346 is a density bounded 
nebula. 



3.1. Constant density models 

The simplest models that can be computed are spherical with constant density (models 
L.l, S.l, and B.l). The density is given by the root mean square density, iV e (rms), and takes 
the value 9.00 cm -3 . This value was derived from the Ha flux from Kennicutt & Hodge 
(1986), a C(Ha) of 0.10 dex (Paper I), a distance to the SMC of 64 kpc (Reid 1999, and 
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references therein), and a radius of 210". 

Using these models we find that most of the predicted line ratios differ considerably 
from the observed ones. In particular the lines of low degree of ionization, such as those of 
Oil and Sn, are very faint (see Table 2). 

By looking at the [O n] and [O ill] line intensities it can be seen that the degree of 
ionization increases in the B.l, A.l, S.l sequence this result permits us to estimate the 
differences in the models due to the ionization field adopted. 

3.2. Models with filling factor different from 1.0 

Considering the filamentary aspect of NGC 346 (e.g. Ye et al. 1991), the densities 
derived in Paper I, and the weakness of the lines with low degree of ionization, we decided 
to compute a model with a filling factor different from 1, and a constant N e (local). 

The filling factor can be estimated by means of the relation: 

A" 2 (rms) = eiV c 2 (local), (2) 

where iVe(rms) is the root mean square electron density, and N c (local) is the electron density 
determined through a forbidden-line ratio or from the He I lines based on the maximum 
likelihood method (see Paper I). The He I lines originate in the (He n) zone, therefore the 
density derived from them is called the (He il) density. 

We computed models of the L series with different filling factors and the best fit to the 
observed O ii/H/3, O m/H/5, and O Ill/O II ratios is provided by model L.4 with an N e (local) 
of 80cm~ 3 and a filling factor e of 0.0127. By comparing the line ratios predicted by model 
L.4 with observations (see Table 2) we find a much better fit than that provided by model 
L.l. Also in Table 2 we present models L.2 and L.3 to show the sensitivity of the fit to the 
adopted density (or filling factor). 

We also computed models of the S series with different filling factors and present them 
in Table 3. Model S.2, which has the same density and filling factor as model L.2, provides a 
considerably better fit to the observations than the S.l model, in agreement with the result 
obtained with the L series. As with the case of the L.l and S.l models, the S.2 model shows 
a higher degree of ionization than the L.2 model. By increasing the density it is possible to 
obtain a better O ll/H/3 fit with a model of the S series, but the fit of the O Ill/O II and the 
Olll/H/5 ratios becomes poorer. 

Models with 0.005 < e < 0.02 provide a much better fit to the observations than models 
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with e = 1.00 (see Tables 2 and 3). 



3.3. Decreasing density models 



We have two determinations of iV e (local), iV e (He ll)= 143 ± 36cm~ 3 and iV e (S ll) = 
50±15cm~ 3 . Since the sulfur emission originates in the outer parts and the helium emission 
all over the H II region, we decided to compute two models with higher density in the center 
and lower in the outer regions. 

Model L.7 has the following decreasing density law: 



where r and N are, respectively, the inner radius and the density at the inner radius. 
We choose a density law with an exponent a=0.315, which gives a value of 180 cm -3 at 
r = 10 19 cm and a value of 70 cm~ 3 at the boundary (r = 2.0 x 10 20 cm). From the 
observational constraint in the Ha flux we obtained that e = 0.013. 

Model B.2 is made of two concentric spherical shells with N e (local) equal to 140 and 
50 cm" 3 respectively and the following radial intervals: 10 16 to 10 20 and 10 20 to 2.0 x 10 20 cm 
respectively. 

Models L.7 and B.2 provide a better fit to the observations than models L.l and B.l, 
but do not provide a significantly better fit than L.4, the main reference model. 



Models L.5 and L.6 were computed with 0.1 dex higher and 0.1 dex lower heavy element 
abundances than model L.4. The differences in the O/H value are about 2a higher for model 
L.5 and 2a lower for model L.6 than the values derived in Paper I. 

Model L.5 produces a modest increase in respect to model L.4 in the 7(3727+5007) / I(Hf3) 
value without reaching the observed value. Alternatively the difference between the predicted 
and the observed 7(4363)//(5007) value becomes larger than in model L.4. 

Model L.6 produces a higher 7(4363)/J(5007) value than model L.4 without reaching 
the observed value. Alternatively the difference between the observed and predicted 7(3727+ 
5007) / I \H (3) value becomes larger than in model L.4. 




(3) 



3.4. Models with different chemical composition 



- 9- 



3.5. Plane-parallel models 

We decided to study a completely different geometric arrangement to see if it was 
possible to get a better agreement between models and observations. In addition to one 
dimensional models with spherical symmetry, it is possible to compute one dimensional 
plane parallel models. For these models we have adopted five observational restrictions 
associated with the geometry: a) the distance to the object, b) the subtended solid angle, 
c) the observed Ha flux, d) an iV e (local)=100 cm" 3 , and e) the incident flux per unit area 
in the direction of the line of sight, that we assume to be constant all over the face of the 
nebula. The incident flux per unit area is given by the ionizing stellar flux divided by the 
area perpendicular to the line of sight. 

The line intensities predicted by one of our plane parallel models (B.3) are presented in 
Table 3. The model has an area perpendicular to the line of sight of 1.26 x 10 41 cm 2 and 
a length in the direction of the line of sight of 0.9 x 10 20 cm, which implies an e of 0.0216. 
The model provides a poorer fit than the spherical models because it shows a considerably 
lower degree of ionization than the observations. 

Other geometries require two or three dimensional photoionization codes and are beyond 
the scope of this paper. 

3.6. Ionization bounded 

To explore further the possibility that NGC 346 is a density bounded Hll region, we 
have studied three line ratios involving lines with low degree of ionization: I (X6300)/ 1 (Hj3), 
I (X6717)/ 1 (H/3), and I (X3727)/ 1 (H/3). These ratios are stronger in the outer zones of an 
ionization bounded H II region than in the outer zones of a density bounded H II region. In 
Table 4 we show these line ratios for each model as well as the observational values from 
Paper I. For I (X6300) / 1 (H (3) , all the density bounded models give a value of the line ratio 
in agreement with the observational limit, and models L.2 to L.6, B.2, and S.2 show good 
agreement with the other two lines. 

In Table 4 we present these three line ratios for two ionization bounded models, IL.4 
and IS. 2. These models are similar to the density bounded models L.4 and S.2 the only 
difference is that the intensity of the ionizing flux is 45% smaller than in the density bounded 
models so that all the ionizing photons are trapped by the nebula. For the IS. 2 model 
the three line ratios are stronger than observed, while for the IL.4 model the predicted 
I (A6300) / 1 (H (3) value is smaller than the upper limit while the other two line ratios are 
stronger than observed. To reach agreement between ionization bounded models and the 
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observations it is necessary to reduce N e (local) (and to increase the filling factor accordingly) 
entering into conflict with the observed iV e (He ll)= 143±36 cm -3 and N e (S ll)= 50±15 cm -3 
values. Of the line ratios presented in Table 4 the best discriminant between the density 
bounded and the ionization bounded models is the I (X6300)/ 1 (H/3) value, therefore it is 
important to obtain a determination of this line ratio and not only an upper limit. The 
previous discussion supports the result that NGC 346 is a density bounded nebula. 

Models for NGC 346 with a covering factor of 0.55 will produce intermediate ionization 
structures to those provided by density bounded and ionization bounded models. From the 
discussion of the previous paragraph it follows that these models will produce better fits to 
the lines of low degree of ionization than the ionization bounded models but poorer fits than 
the density bounded models. 



4. Ionization Correction Factors and Abundances 

For most of the elements we observe only one or two stages of ionization, therefore to 
obtain the total abundances relative to H we need to determine the fraction of a given element 
present in the unobserved stages of ionization. Therefore to obtain the total abundances we 
have made use of the following equations: 



N(0) _ N(0 + ) + iV(0++) 
JV(H) ~ iV(H+) ' 

iV(N) 



(4) 



= JCT(N) iV(N+)/iV(H + ), (5) 
= JCF(Ne) iV(Ne ++ )/iV(H + ), (6) 
JCF(S) iV(S + + S ++ )/iV(H + ), (7) 



JV(H) 
iV(Ne) 

W) : 

7V(H) 

im - ICF{kI) — irw) — ' { ) 

^ = JCF(Fe) N(Fe ++ )/N(R + ); (9) 

where the ionization correction factors, ICF's, were estimated from the models presented in 
this paper, see Table 5, and the ionic abundance ratios were taken from Paper I. To obtain 
the N(Fe ++ )/N(H + ) value we made use of the observed intensity of A 4568 in region A and 
its predicted intensity by the models. The total abundances are presented in Table 6. 



The He/H value is given by (e.g. Sauer & Jedamzik 2001; Peimbert 2001 and references 
therein): 



iV(He) / N e N(tie°)dV + / N e N(Re + )dV + J N e N(Re ++ )dV 
N(R) ~ J N e N(R°)dV + J N e N(R+)dV ' 

= ICF{Re) fN e N(H+)dV ' (10) 

where ICF(Re) is the helium ionization correction factor. To determine a very accurate 
He/H ratio a very precise value of ICF(He) is needed. As expected for our density bounded 
models it is found that ICF(He) = 1.000, because there is no transition zone from ionized 
to neutral He and H regions. 



5. Discussion and Conclusions 

We have presented a set of photoionization models for NGC 346. The advantage of these 
models relative to those models present in the literature of giant extragalactic H II regions 
is that we know the spectra of most of the ionizing stars and the absolute V magnitude for 
all of them, therefore we did not have to assume the IMF, the age of the cluster, the star 
formation rate, nor the fraction of ionizing photons that are trapped by the nebula. 

The models were computed with different geometries, different density distributions, 
different chemical compositions and different radiation fields. We did not compute models 
composed via a superposition of separate individual Stromgren spheres because a consid- 
erable amount of the ionizing flux escaped from the nebula implying that the Stromgren 
spheres produced by each star overlap, and that the single Hll region model is a good 
approximation. 

By comparing the models with the observed line intensity ratios we obtain the following 
results: 

1) Spherical models with constant density, with e = 1.00 (models L.l, S.l, and B.l), 
attain a high degree of ionization that is not observed. Moreover, these models have an 
iV e (rms) = 9.00 cm" 3 , this density is smaller than those derived from the [Sll] and He I lines 
(see Paper I). 

2) Spherical models (L.4, L.7, S.2 and B.2) can reproduce the observed degree of ion- 
ization and show good agreement with the [S il] and He II densities. 
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3) Models L.7 and B.2 include more realistic density distributions, but their predictions 
are not very different to those of models L.4 and S.2. 

4) The /(Hell, 4686)/J(H/3) ratio is very small and is very sensitive to small changes in 
the radiation field and the geometry of the Hn region. Model L.l has the same geometry 
as models S.l and B.l but shows a difference of about an order of magnitude in the /(Hell, 
4686)//(H/3) ratio. The differences are due to the three radiation fields used which vary in 
the fraction of photons with an energy higher than four Rydbergs, while in model B.l we 
use a blackbody spectrum as an ionizing source, in models L.l and S.l we use two different 
sets of stellar atmopheres (see Figure 1). 

Models L.2 to L.6 have similar /(He II, 4686) / /(H/3) ratios while model L.7 shows a ratio 
that is about two orders of magnitude fainter than the other L models. All the models of 
the L series have the same radiation field, in this case the difference is due to the geometry, 
the L.7 model has an inner hole of 10 19 cm while the other L models have an inner hole of 
10 16 cm, the different geometry dilutes the radiation field in the inner regions of the nebula 
causing the drop in the helium ionization degree. A similar drop in the ionization degree, 
due to geometrical effects, is present in the plane parallel model B.3, which shows an /(Hell, 
4686)//(H/3) ratio two orders of magnitude fainter than the spherical models B.l and B.2. 

5) Plane parallel models present a considerably lower degree of ionization than is ob- 
served in NGC 346. For example they underestimate the /([Ariv], 4740) //([Arm], 7135) 
and the /([Olll], 5007) //([On], 3727) ratios. On the other hand spherical models (like L.4, 
L.7, S.2 and B.2) present line ratios in good agreement with the observations. 

6) None of the models are able to reproduce the observed /(3727 + 5007)//(H/9) ratio 
nor the /(5007)//(4363) ratio. This is a crucial result which implies that photoionization 
models predict lower temperatures than those observed. This result probably implies that 
there are additional heating sources not considered by the photoionization models. Similar 
results were obtained by Stasihska & Schaerer (1999) for I Zw 18, Luridiana, Peimbert, & 
Leitherer (1999) for NGC 2363, and Luridiana & Peimbert (2001) for NGC 5461. 

7) The total abundances of N, Ne, S, Ar, and Fe were obtained from the ICF's predicted 
by model L.4. By comparing the empirical ICF's in Paper I with those derived from 
model L.4 we find excellent agreement for Ne and Ar, good agreement for S, and only a fair 
agreement for N. Based on empirical methods and on photoionization models it has been 
found that the + zone coincides with the N + zone and consequently that the N/O ratio is 
equal to the N + /0 + ratio (e. g. Peimbert & Costero 1969; Garnett 1990; Mathis & Rosa 
1991; Stasinska & Schaerer 1997, Paper I); alternatively the N abundance derived from the 
NGC 346 CLOUDY models differs by about a factor of two from that derived from the N+/0+ 
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ratio (see Table 6), Stasiriska & Schaerer (1997) have studied this discrepancy and conclude 
that it is due to the different ionizing spectra used by different authors. 

From a detailed comparison between stellar ionizing photon input and the observed 
Ha flux we infer that some 45% of the photons produced by the ionizing stars escape from 
NGC 346. This result is independent of the geometry and of the density distribution of 
the nebula, and implies that this object must be a major source of ionizing flux for the 
surrounding diffuse interstellar medium. The ionization bounded models we have tried fail 
completely to match the intensities of lines with a low degree of ionization, but these can 
be matched using density bounded models, which supports strongly the conclusion that 
NGC 346 is density bounded. 
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Table 1. Spectral Types, Temperatures, Mass, and Ionizing Fluxes 



Spectral Type 


T 


M 


Number of stars 


Ionizinj 


y Flux 




(K) 


(M©) 






(10 39 erg s" 1 ) 


03V 


51230 a 


51.3 a 


1 




2.02 




04V 


48670 a 


44.2 a 


1 




1.43 




05.5V 


44840 a 


35.5 a 


3 


(1) 


2.46 


(0.82) 


kJO V 


4O00U 


OO. 1 


3 


(1) 


2.07 


(0.69) 


06.5V 


42280 a 


30.8 a 


2 


(1) 


1.14 


(0.57) 


07V 


41010 a 


28.8 a 


3 


(2) 


1.40 


(0.94) 


07.5V 


39730 a 


26.9 a 


2 


(2) 


0.76 


(0.76) 


08V 


38450 a 


25. l a 


6 


(7) 


1.88 


(2.20) 


08.5V 


37170 a 


23.6 a 


1 


(11) 


0.26 


(2.82) 


09V 


35900 a 


22. l a 


2 




0.43 




09.5V 


34620 a 


20.8 a 


5 




0.86 




05.51 


43210 a 


45.4 a 


1 




4.93 




071 


38720 a 


37.4 a 


1 




1.77 




WN4(5980A) 


52000 b 


18.0 C 


1 




7.01 




O7Ia(5980B) 


31135 b 


37.4 a 


1 




2.86 




Total 






33 


(25) 


31.28 


(8.80) 



a Vacca et al. (1996). 
b Schweickhardt & Schmutz(1999). 
c Koenigsberger et al. (1998) 



Table 2. Photoionization models based on a set of atmospheres with = 0.2Z Q 



Line ratios a 


Reg. A 


Reg 3 and 13 


Mod L.l 


Mod L.2 


Mod L.3 


Mod L.4 


Mod L.5 


Mod L.6 


Mod L.7 


[Om] J(A5007)/[On] /(A3727) 


0.738 


0.750 


1.633 


0.492 


0.383 


0.595 


0.597 


0.578 


0.769 


[Sin] 7(A6312)/[Sn] /(A6725) 


-0.857 


-0.744 


-0.004 


-0.883 


-0.974 


-0.802 


-0.847 


-0.776 


-0.659 


[Anv] 7(A4740)/[Ann] 7(A7135) 


-1.111 


-1.392 


-0.614 


-1.233 


-1.292 


-1.179 


-1.214 


-1.155 


-1.136 


[Om] 7(A4363)/i"(A5007) 


-1.854 


-1.915 


-2.041 


-2.015 


-2.010 


-2.021 


-2.100 


-1.951 


-2.018 


He II /(A4686)/Hei /(A4471) 


-1.155 


-1.430 


-0.905 


-1.286 


-1.337 


-1.242 


-1.243 


-1.243 


-3.209 


[Sir] 7(A6716)/7(A6731) 


0.144 


0.142 


0.152 


0.116 


0.106 


0.124 


0.124 


0.123 


0.123 


[S n] /(A4069 + 4076)//(A6725) 


-1.115 


-0.940 


-1.104 


-1.050 


-1.038 


-1.060 


-1.086 


-1.037 


-1.060 


[Om] I (X5007)/ 1 (H/3) 


0.735 


0.710 


0.625 


0.541 


0.518 


0.559 


0.579 


0.526 


0.609 


fOnl I(XX72T\II(Hff\ 


-0.003 


-0.040 


-1.008 


0.049 


0.135 


-0.036 


-0.018 


-0.052 


-0.160 


He II I(\4686)/I(Hf3) 


-2.571 


-2.844 


-2.419 


-2.804 


-2.856 


-2.759 


-2.753 


-2.766 


-4.708 


L(Ha)(dcx) 






38.987 


38.978 


38.987 


38.977 


38.988 


38.981 


38.989 


-/V e (rms)(cm~ 3 ) 






9.00 


9.00 


9.00 


9.00 


9.00 


9.00 


9.00 


radius (10 20 cm) 






2.00 


2.01 


2.01 


2.01 


2.00 


2.04 


2.02 


N e (local) 






9.00 


100 


130 


80 


80 


80 




€ 






1.00 


0.0081 


0.005 


0.0127 


0.0127 


0.0127 


0.013 


[o] b 






8.11 


8.11 


8.11 


8.11 


8.21 


8.01 


8.11 


T e (rad)K 






11200 


11400 


11500 


11400 


10700 


12000 


11400 


T e (vol)K 






10800 


11200 


11300 


11100 


10600 


11700 


11100 


T e ([OlIl])K 


13070 


12430 


10800 


11300 


11300 


11200 


10600 


11800 


11200 



a Given by log J(Ai)/J(A 2 ). 

b Gaseous abundances given by 12 + log A^OV-N^H). 



Table 3. Photoionization models based on other sets of stellar atmospheres 



Line ratios a 


Reg. A I 


leg 3 and 13 


Mod S.l b 


Mod S.2 b 


Mod B.l c 


Mod B.2 C 


Mod B.3 C 


[Om] 7(A5007)/[On] /(A3727) 


0.738 


0.750 


1.918 


0.716 


1.435 


0.759 


0.190 


[Sin] J(A6312)/[Sn] /(A6725) 


-0.857 


-0.744 


-0.065 


-0.927 


-0.054 


-0.587 


-0.882 


[Anv] 7(A4740)/[Arin] /(A7135) 


-1.111 


-1.392 


-0.346 


-0.994 


-0.812 


-1.156 


-2.581 


[n tttI t / \ aqrq\ I t ( \ K(\r\ r 7\ 
[U IIIJ 1 {A4obo)/l (AoUUf ) 


-1.854 


-1.915 


-1.981 


-1.976 


-2.077 


-2.039 


-2.078 


He II 1 (A4Doo)/rlc I i(A4471) 


-1.155 


-1.430 


0.293 


0.183 


0.359 


0.359 


-1.916 


[So] J(A6716)/i (A6731) 


0.144 


0.142 


0.150 


0.117 


0.152 


0.112 


0.116 


[S n] /(A4069 + 4076)//(A6725) 


-1.115 


-0.940 


-1.080 


-1.043 


-1.125 


-1.065 


-1.065 


[O in] / (A5007)/i (-Hp) 


0.735 


0.710 


0.663 


0.607 


0.595 


0.585 


0.478 


[Oil] l(M727)/l(Hp) 


-0.003 


-0.040 


-1.255 


-0.109 


-0.840 


-0.174 


0.288 


Hen /(A4686)//(,H'/3) 


-2.571 


-2.844 


-1.255 


-1.360 


-1.132 


-1.137 


-3.338 


L(Ha)(dcx) 






38.981 


38.973 


38.975 


38.989 


38.993 


A^e (rms) ( cm ~~ 3 ) 






9.00 


9.00 


9.00 


9.00 


9.00 


radius (10 20 cm) 






2.01 


2.01 


1.94 


2.00 




R x (10 20 cm) 












1.01 




N e (local) 






9.00 


100 


9.00 


(140,50) 


100 


e 






1.00 


0.0081 


1.00 


0.0175 


0.0216 


[o] d 






8.11 


8.11 


8.11 


8.11 


8.11 


T e (rad)K 






11700 


11700 


10800 


11100 


10900 


T e (vol)K 






11200 


11400 


10400 


10800 


10900 


T e ([OlIl])K 


13070 


12430 


11300 


11500 


10400 


10800 


11000 



a In units of log /(Ai)/7(A 2 ). 
b Z* = Z Q 

c Blackbody ionizing spectrum. 

d Gaseous abundances given by 12 + log N(0)/N(H). 
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Table 4. Line ratios of [0 i] , [S il] , and [0 il] 



Models 7(A6300)//(#/3) 7(A6717)/7(///3) 7(A3727)//(///3) 

= O.2Z 

L.l -5.402 -2.244 -1.008 

L.2 -3.155 -1.140 0.049 

L.3 -2.953 -1.042 0.135 

L.4 -3.340 -1.230 -0.036 

L.5 -3.289 -1.279 -0.018 

L.6 -3.362 -1.174 -0.052 

L.7 -3.609 -1.327 -0.334 
Blackbody Spectrum 

B.l -5.277 -2.136 -0.840 

B.2 -3.878 -1.433 -0.174 

B.3 -3.328 -1.029 0.288 
Z*. = Zq 

5.1 -5.600 -2.297 -1.255 

5.2 -3.219 -1.116 -0.109 
IL.4 a -2.270 -0.861 0.234 
IS.2 a -1.707 -0.844 0.193 

Observations 13 < -2.016 -1.284 -0.040 



ionization bounded models, see text 
b For Region A and from Paper I 
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Table 5. Ionization Correction Factors for Region A 



Element 


Paper I 


This paper 


N 


5.62 


11.22 


Ne 


1.23 


1.29 


S 


1.98 


1.29 


Ar 


1.21 


1.22 


Fe 




8.00: 
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Table 6. Total Abundances a 



Element 


NGC 346 


NGC 346 


Sun b 


Orion c 


M17 d 




Paper I 


This paper 








N 


6.51 


6.81 


7.92 


7.78 


7.90 


C 


8.15 


8.15 


8.83 


8.72 


8.87 


Ne 


7.30 


7.32 


8.08 


7.89 


8.02 


S 


6.59 


6.40 


7.33 


7.17 


7.31 


Ar 


5.82 


5.82 


6.40 


6.49 


6.60 


Fe f 




6.16 


7.50 


6.11 


6.69 



a In units of 12 + Log N(X)/N(K). 
b Grevesse & Sauval (1998). 
c Esteban et al. (1998). 

d Peimbert, Torres- Peimbert, & Ruiz (1992); Esteban et al. 
(1999). 

e Gas plus dust for the Hn regions. 
f Gaseous content only. 
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1 2 3 4 5 

Energy(Rydberg) 



Fig. 1. — Ionizing continuum spectra used in the models. All the spectra are normalized to 
the same total number of ionizing photons. 




Fig. 2. — IMF derived from Table 1. The straight line represents the IMF derived from 
the four central bins and it has a slope of -3.7±0.4. This IMF has an upper mass cutoff of 
54.1 M and it is only representative for masses higher than 24 M . 



